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1. Introduction 



The properties of a relativistic particle in a one-dimensional periodic structure play an 
important role in understanding many phenomena in Solid State and Nuclear Physics 
For this purpose the relativistic generalization of the classic Kronig-Penney 
model mis], the numerical solution of the Dirac equation P and a class of exactly 
solvable nonlocal separable periodic potentials 0] have been studied. Nevertheless, 
due to its elementary character, the Dirac Kronig-Penney model remains the most 
frequently used model for studying one-dimensional relativistic band structures. 

A new method for constructing exactly solvable non-relativistic periodic 
potentials has been recently proposed [S]. Starting from an easily solved (non- 
periodic) potential, a new solvable potential is generated by using the familiar Darboux 
transformation. This new potential is then considered on a bounded interval and 
its periodic extension is constructed. In this paper we generalize this procedure to 
construct a solvable scalar potential for the Dirac equation. 



2. Band structure for a periodic potential 

Here we recall the main features of the Dirac equation with a periodic potential (see 
e.g. 011). 

Consider the one-dimensional Dirac equation 

h%l;{x) ^ Ei:(x), h = iayda, + V{x), (1) 

where V{x) is a real self-adjoint matrix potential, which we assume has the canonical 
representation (for more details see e.g. |S]) 

V{x) = p{x)az + q{x)ax ■ (2) 

'^x.y,z are the usual Pauli matrices and ipi^) — {ipi{x),'ip2{x)y (the superscript "f 
meaning the transpose) is a two-component spinor wave function. For any real value 
of the parameter E equation Q has two real linearly independent solutions ipi^) and 
(p{x). If E is chosen to be real, then the wave functions can be chosen real as well. 
The adjoint equation for 1^9(2;) is 

- id^Lp+{x)ay + ip+{x)V{x) = Eip+{x) . (3) 

If we multiply by the row- vector (x) from the left and JSJ by the column- vector 
ipix) from the right, the difference takes the form [H^((/3, ^)]a; — (the subscript x 
stands for the derivative with respect to x) meaning that the function 

W{(p, ijj) = (p+ {x)icryip{x) . (4) 

is (a non-zero) constant. Hence, W{(p,ip) here plays the role of a Wronskian for the 
Dirac equation. Since ip and ip form a basis set for the two-dimensional solution space 
of for any given E, any solution of with the same E may be expressed as 

Xix) = A^^ix) + Bifix) , (5) 

where A, B € C Without loss of generality we may assume that ip and ip have been 
chosen so W{ip, ip) = 1. 

Suppose now that V{x) is a periodic potential: 

V{x) ^Vix + T) TeK. (6) 
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By Floquet's theorem, in this case has Block solutions such that 

x{xo)^ f3x{xo + T) /3 = const eC, (7) 

where xq is any point on the real axis. By using |(SJ we get the following equation for 
the coefficients A and B: 

^{xo)A + ^{xo)B = P[^/j{xo + T)A + ^{x^ + T)B] . (8) 

which is equivalent to the system of linear homogeneous equations: 

[V-iCxo) - /3^i(a;o + T)]A + [^i(a;o) - /3(^i(a;o + T))]B = 

[iP2{xo) - PMxo + T)]A + [ip2{xo) - f3Mxo + T)]B = ' 

For such a system to have a non-trivial solution the determinant of coefficients must 
vanish, which is the case only if (3 satisfies the quadratic equation 

/J^ + 1 = 0, (10) 

where 

D = W{i^{xa + T), v?(.To)) + W^(V'(xo), (^(xo + T)) (11) 

and we have used the property W{'4>,if) — 1. We get from (|10|) the two possible values 
for (3 

/3i,2 = i?/2± V^V4-1, /3i/32 = l. (12) 

This result is quite similar to that for the Schrodinger equation with a periodic 
potential The function D ~ D{E) is called the Lyapunov function or Hill 

determinant for the Dirac equation, and is real since ip and ip are. 

It is clear that for |Z)(£')| < 2 the values of /3 lie on the unit circle of the complex 
/3-plane, so we can write them 

Pi = exp {iKT) , /?2 = exp {~iKT) . (13) 

with real K. It follows from Q that the Bloch solutions are bounded on the whole 
real axis. Hence, all real E with \D{E)\ < 2 lie in allowed bands. Since for |D(£^)| > 2, 
/3i_2 are real and positive, either Bloch solution is unbounded. All such E form the 
forbidden bands, while the band edges are given by |Z?(i?)| — 2 and in general form 
the unbounded sequence 

. . . E^4 < E^3 < E^2 < E^i < El < E2 < Es < E4.. . . (14) 

An important feature of the band structure for the Dirac equation is that allowed 
bands may exist for negative energy. 



3. Darboux transformation operator 

The Darboux transformation provides a very powerful method for finding new exactly 
solvable potentials, both for the Schrodinger equation j7j and the Dirac equation 
In this section we shall describe the main ideas for this method following (S) |5| . 

Consider the stationary one- dimensional Dirac equation for the Hamiltonian Iiq: 

ho^j{x) = {layd, + Vo{x))^{x) = E^{x) . (15) 

Suppose the solutions of (|15|l arc known for all E and we wish to solve a second Dirac 
equation having Hamiltonian hi. 

hiip{x) = [iaydx + Vi{x))ip{x) = Eip{x) . (16) 
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Instead of solving H16|l directly, one looks for a transformation operator (or intertwiner) 
L such that 

Lho = hiL. (17) 

If such an operator can be found, then the eigenspinors of hi may be obtained by 
applying L to those of ho: ip{x) = Lip{x). 

It has been shown (see [HI El for details) that a first order differential operator 

L ^ ~ u {x)u-^ (x) (18) 
satisfies (|17f) with 

Viix) = Va + [lay, u'ix)u-\x)], (19) 
where u{x) = {u'^^\x) , u^^^x)) is a 2 x 2 matrix satisfying the equation 

hou{x) = u{x)A, A = diag(Ai, A2) . (20) 

It is easy to check that the spinors u^^\x), j = 1,2 are eigenfunctions of the 
Hamiltonian /iq corresponding to the eigenvalues A^. The matrix u{x) is called a 
transformation function. 

If -E 7^ Ai or A2, then by applying L to any solution to H15|l . we obtain a 
corresponding solution to H16|l : if £' = Ai or A2, then Lu^^'^\x) — 0. Nevertheless, in 
the latter case it can be shown that the function 

vix) = iu\x)r' (21) 

is a matrix eigenfunction for the Hamiltonian hi with eigenvalue matrix A: hiv{x) ~ 
v(x)h^. Therefore the spinors v^-^'>(x\ j = 1, 2 are eigenfunctions of hi with eigenvalues 
Aj. Solutions to (16), v^^''^^{x) for E = \i and E = A2, may be found by means of the 
properties 

W{v'^^\x),v'^^\x)) = \, j = l,2. (22) 

Therefore, with the help of the Darboux transformation one is able to obtain the 
solutions to the transformed equation \V6\ for any value of E provided the solutions 
to the original equation H15|) are known. 

4. Darboux transformation for a scalar potential 

A scalar potential is specified by a single function pq{x) = m + S(){x), x G R, so 

Vo = im + Soix)).. = ( ™ ) , (23) 

where m is the particle mass. For our purposes, it is more convenient to use the 
alternative form for the scalar potential 

/ , r. ^ f m + 5*0(2;) \ 

Vo = {m + So{x))a^=^^^_^g^^^^ j' ^^4) 

These two potentials are related by the unitary transformation 

Vo^U-'VoU,^j^U4', (25) 

where 

U ^{l + iay)/V2. (26) 
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In general, when a Darboux transformation is applied to a scalar potential, the 
resulting potential is not scalar, but below we formulate additional conditions that 
will prevent this (see [S] for details) . 

If the spinor u^^^ = (itii,U2i)* is a solution to the Dirac equation (|15|) for the 
potential H24() . for some energy E — \,it is easy to check that the spinor u'^^ = — cr^u'^^ 
is a solution to the same equation for energy E = —A. The transformation matrix 
u{x) constructed from the spinors u^^^ and w'^-* can be seen to give 

= ( „). (27) 



(lnu2i)' 
and the transformed potential H19|l takes the form 

= {m + Si{x))a, , (28) 
Si{x) = So + {lnu2iy - {InunY (29) 

which is clearly scalar. 

With the help of H18() and H27() one finds the solutions to the transformed equation: 

T I ( ■'/''i - (lnuii)'-0i \ , , 

Lp^L-^^\ ; \ ' \. (30) 
\W2- (lnu2i) -02 ) 

Note that the components of this spinor are given by the same expressions that 
appear in the Darboux transformation, with transformation functions mh, U21, for 
the Schrodinger equation (see e.g. [TT]). 

5. New periodic scalar potential 

The free particle Dirac Hamiltonian can be written in terms of the scalar potential 
H24|l with ^0 = 0. One can easily verify that the spinors 



u^^^(x) = (cosh(7x — a) , cosh(7x + a))*, e^" — ^ — — — , 

V "^ + 7 (31) 

u'^^^x) = -a,u'-^\x) 



are eigenspinors of ho corresponding to the eigenvalues A1.2 = ±A = ±\/m^ — 7^. By 
constructing the transformation function from these spinors and using (|29|l , we arrive 
at the reflectionless (one-soliton) potential corresponding to 

27^ 

= ^x\o . (32) 

TO + A cosh 27a; 

which was found previously |12) . It is easily checked that the two spinors given by 
the matrix (u+)~^ are square integrable. Hence the potential H28|) . (|32|) has the two 
discrete levels E — ±X. 

Linearly independent solutions of the Dirac equation for potential H32|l found by 
means of ()30|l are 

7/ . _ E / coskx - lwi{x)smkx \ , . 

\ cofi{kx — 6) — ^W2{x) sm{kx — 6) y ' 

, _ 1 / kainkx + wi{x) coskx \ , , 

^ ' fc2 ksm(kx — S) + W2{x) cos{kx — S) y ' ^ ' 

where k = V E'^ — m^, 6 = arctan(fc/TO), wi^2ix) = 7tanh(7a; ^ a). These are 
normalized such that their Wronskian is unity: W {^e{x) , ifiEix)) = 1. 
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Figure 1: Periodically continued scalar potential. 



Let us now restrict the potential as given by (|32|l to the interval [—a, a], a > 
and continue it to the entire a;-axis by periodicity to obtain the continuous potential, 
having period T = 2a, shown in Fig.l for m — 2, \ — \^ and a — 1. Having foimd the 
solution to the Dirac equation allows us to calculate the Lyapunov function Hll|l . For 
this purpose we first have to calculate Wronskians at xq = —a. So, using H33|l and 
H34|l after some simple algebra we find for the first one the expression: 



W(^(a),^(-a)) 



E 



+k cos{2ka) sm{5) — W2{a) ( cos(2fca — S) 



wii^a) ( cos(2fca + S) 
1 



■ cos(fca + S) sin(fca) i«i(a) 



cos(fca) sin(fca — 6) W2{a] 



(35) 



In the derivation we have used the symmetry relations wi(— a) = W2(a), W2{—a) = 
Wi (a) . After calculating the second Wronskian in (|f 1|) we obtain 

D{E) = t: \2wi{a) cos(2fca + 5)- 2w2(a) cos(2fca - S) 

+ 7^ '■ 

+ i(fc2 - wl(a)) sin(2fca + S) - ^(k^ - w?(a)) sin(2/ca - 5)] 
k k 

The behavior of this function is sketched in Fig. 2 



(36) 









-2 



Figure 2: Lyapunov function for a periodic scalar potential. 

for TO = 2, A = 1, and a — 1 and leads to the following values for the lowest band 
edges in the positive spectrum: Eq = 0.738, Ei = 1.381, E2 = 2.164, E3 = 3.274, 
Ei = 3.335, E5 = 4.802, Eq = 4.827, E7 = 6.352. Since D{E) is an even function of 
the energy, the negatives of these values form the band edges in the negative spectrum. 
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We conclude this section by indicating how the actual energy bands may be 
constructed from our results. From {Tj), and since the Wronskian has the 

value unity, we see that 

cos(2i4:a) = ^D{E) (37) 

so that the wave vector K is obtained as an explicit function of the energy in the j-th 
energy band by 

K = ^ arccos ^D{E) for E-j^i < E < Ej . (38) 

By symmetry it is sufficient to consider only positive as we have implicitly done 
here. In general the inversion of this to get i? as a function of K (the so called 
dispersion law) must be carried out numerically, but is trivially done graphically. The 
lowest positive energy band for the model above is shown in Figure 3. 




0.2 0.4 0.6 



Figure 3: Lowest energy band for the case 7 = 1, m = 2 for tlie scalar potential given by 1321 . 



6. Discussion 

In this note, using a specific example of one-soliton potential, we have shown how 
a great variety of exactly solvable continuous periodic local potentials for the one 
dimensional Dirac equation may be constructed. Furthermore, this procedure may 
be manipulated so that only elementary mathematical functions are involved in the 
construction and in the solution. Up till now, the only examples available have been 
simple variants of the classic Kronig-Penney model, yet this has been a key test bed 
in the relativistic theory of surface states ^2!) for example. We feel that the class of 
models presented here will be equally useful and lead to new insights in this and other 
areas. 
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